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Abstract 

We study an upper bound on masses of additional scalar bosons from the electroweak precision data and 
theoretical constraints such as perturbative unitarity and vacuum stability in the two Higgs doublet model 
taking account of recent Higgs boson search results. If the mass of the Standard-Model-like Higgs boson is 
rather heavy and is outside the allowed region by the electroweak precision data, such a discrepancy should 
be compensated by contributions from the additional scalar bosons. We show the upper bound on masses 
of the additional scalar bosons to be about 2 (1) TeV for the mass of the Standard-Model-like Higgs boson 
to be 240 (500) GeV. 
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The Standard Model (SM) for elementary particles based on the SU{3)c x SU{2)i x U{1)y 
gauge group has been tested accurately |jl|. However, no Higgs boson has been confirmed yet. 
Discovery of the Higgs boson is the most important issue at the Fermilab Tevatron and the CERN 
Large Hadron Collider (LHC). Direct searches for the Higgs boson at CERN LEP have set a lower 
mass bound on the SM Higgs boson to be 114.4 GeV[^. The Tevatron experiment has excluded 
the mass of the SM Higgs boson around 160 GeV[^. Recently, the first results from the ATLAS 
and CMS experiments at the LHC have been reported |, |. The Higgs boson mass around 160 
GeV and 300-450 GeV has been excluded by the data with the integrated luminosity of about 1 
fb-i. 

It is well known that an upper bound on the mass of the Higgs boson is obtained by the tree 
level unitarity for elastic scattering processes of longitudinally-polarized vector bosons, such as 
W^W£ —7- W^W£ . In the SM, since the scattering amplitudes are proportional to the Higgs 
boson mass in the high energy limit, a large Higgs boson mass leads to a violation of the unitarity. 
Consequently the upper bound is obtained on the mass as about 710 GeV(|, 0. On the other hand, 
if the Higgs boson is absent, the scattering amplitudes grow for high energies. The violation of the 
tree level unitarity then occurs at y/s ~ 1.2 TeV, where ^/s is the centre-of-mass energy of the WW 
scattering[^. The LHC Higgs search experiment is expected to cover the entire range of the SM 
Higgs boson mass. Even if the Higgs boson is not found, some new physics beyond the SM must 
show up below the TeV scale. If we introduce the cutoff scale A into the model, more sensitive 
upper and lower bounds are obtained on the SM Higgs boson mass as a function of A[|-|ll|. 

From the electroweak precision data with the theoretical study for radiative corrections [^], the 
mass of the Higgs boson in the SM is indicated to be rrih = 90j^22 GeV and nih < 161 GeV at the 
95% Confidence Level (CL)|l^]. Notice that this indirect bound on the mass cannot be applied if 
new physics exists below the TeV scale and affects the calculation for the radiative correction. In 
such a case, even if the Higgs boson is heavy, the electroweak precision data can be satisfied by 
the contribution from the new physics. 

In the two Higgs doublet model (THDM), radiative corrections to the electroweak observables 
have already been calculated, and the possible allowed regions for the parameter space are evaluated 
under the electroweak precision data[14, ^ and theoretical constraints! 16 -2C]. Flavor physics data 



such as 6 — )• S7[^, 22], i? — )• Tz^[g^] and tau leptonic decays] 24, ^ in the THDM can further 



constrain the parameter space depending on types of Yukawa interactions. In particular, the mass 
of charged Higgs bosons is bounded from the 6 — )• 57 data to be greater than 295 GeV[^ by 
assuming the Type-II Yukawa interaction. 
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In this Letter, in light of recent Higgs boson searches, we reanalyze the constraint on the 
parameters in the THDM by using the electroweak precision data and the theoretical constraints 
from tree level unitarity and vacuum stability. In particular, we show an upper bound on the 
masses of the additional heavy scalar bosons can be obtained in the THDM depending on the mass 
of the SM-like Higgs boson. For a relatively large mass of the SM-like Higgs boson, a relatively 
large mass difference between the CP-odd Higgs boson and the charged Higgs boson is required in 
order to satisfy bounds from the electroweak precision data. They are bounded from above by the 
theoretical constraints. For a SM Higgs boson mass to be 240 (500) GeV, an upper bound on the 
mass of the CP-odd Higgs boson of about 2 (1) TeV is obtained. 

The most general THDM with the doublet fields <I>i and $2 are constrained by flavor changing 
neutral current (FCNC) processes. We here consider the model with the softly-broken discrete Z2 



symmetry under <I>i — )■ <I>i and ^2 — ^ —^2 to avoid FCNC constraints] 27]. There are four kind of 



Yukawa interaction under the discrete symmetry|21]. In this Letter, we do not specify the type of 



the Yukawa interaction because it does not affect the following discussions. The Higgs potential is 
then given by 

VTHDM ^ +^2^t$^ + ^2^t^2 _ ^2 j^^J^^ ^ ^t^^^ ^ + y (^^2^2)' 

+ A3(^>I$l)($5$2) + A4($l$2)(^^2^l) + y [(^1^2)' + {^\^lf\ ■ (1) 

The soft-breaking mass parameter mg and the coupling constant A5 are complex in general. We 
here take them to be real assuming that CP is conserved in the Higgs sector. 

The Higgs doublets = 1,2) can be written in terms of the component fields as 



where the vacuum expectation values (VEVs) v\ and V2 satisfy \A^^~+~*¥ = v 246 GeV. The 
mass eigenstates are obtained by rotating the component fields as 





(3) 



Ml \h I \z2j \AJ \UJ2 

where uj^ and z are the Nambu-Goldstone bosons, h, H, A and are respectively two CP-even, 
one CP-odd and charged Higgs bosons, and 

/ cos y — sin ( , 

m ={ I • (4) 

\ sin 9 cos ( 
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The eight parameters m\-m^ and A1-A5 are replaced by the VEV the mixing angles a and 
/?(= tan^-*^ the Higgs boson masses m^, mjj,m^ and m^^, and the soft Z2 breaking parameter 
= m3/(cos /3 sin/3). In particular, the quartic coupling constants are expressed in terms of 

physical Higgs boson masses, mixing angles and the soft Z2 breaking mass parameter as 

1 

v"^ cos^ P 



—M sin P + rriTT cos a + m^sin a), (5) 

1 



A2 = ^5 5 — (— cos^ /? + sin^ a + m\ cos^ a) , (6) 

sin p 



A3 



.M2 + (m|,-m2)||i0+2m2,+ 



(7) 



^4 = ^(^^ + "iA-2m^+), (8) 
1 



A5 = ^(M2-mi). (9) 

The coupling constants of the CP-even Higgs bosons with the weak boson hWW and HWW are 
proportional to sin(/3 — a) and cos(/3 — a). When sin(/3 — a) = 1, only /i couples to the gauge 
bosons and behaves as the SM Higgs boson. We call this limit as the SM-like limit. 

As discussed in Ref. [28|, the masses of the heavy Higgs bosons {H, A and H^) are expressed 
for M > t; by 

m| ~ m2 + XiV^[+OiXv^/M% (10) 

while the mass of h is the SM-like form ~ AjV^. When ^ XiV^ the heavier Higgs bosons 
have the common mass ~ M. In this case, the effect of these bosons decouples in the large mass 
limit and the low energy theory becomes the SM with h being at the electroweak scale as the SM 
Higgs boson. On the contrary, when ~ AjW^ the effect of these bosons does not decouple, 
and so-called nondecoupling effects appear in the low energy observables. Notice that the mass 
difference between the heavy Higgs bosons is independent of M, so that the effect on the low 
energy observables can be large if the mass differences are not small. We also note that the 
mass difference between heavy Higgs bosons are related to the violation of the custodial SU{2) 
symmetry [29|], which causes significant deviation in the electroweak rho parameter from the SM 
prediction in the positive direction. As we see soon below, this positive contribution to the rho 
parameter (or the T parameter |jl^) can be used to compensate the negative contribution of the 
heavy SM-like Higgs boson. 

New physics effects on the electroweak oblique parameters are parameterized by the S, T and 
U parameters [12]. By fixing U = 0, the central values of S and T are given by|j|] 

S = 0.03 ± 0.09, T = 0.07 ± 0.08, {pgr = 0.87), (11) 
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where pgrp is the correlation parameters for the analysis. The origin 5 = T = corresponds to 
the SM prediction for the reference value m/j = 117 GeV. The other SM parameters are chosen as 
= 0.23124 ± 0.00016, = 0.01183 ± 0.0016, mt = 173 ± 1.3 GeV and Gp = 1.16639 x 10^^ 
GeV-2. 

In the THDM, the contributions to the electroweak parameters from the scalar boson loops are 



given by 1 14 1 



5*$ = - -^[F Aim H±,mjj±) - sin2(/3 - a)F^{mH,mj^) - cos^(/3 - Q)F^(m/j, m^)] , (12) 



iovr ^ 

+ cos^(/3-a)[FA(m^,m^) - FA(m;„ |, (13) 

where 

Fa mo, mi = Fa mi, mo = \ ^ - ,° % ln-^, 14 

2 mg — m| mf 

X r/1 ^2i^ ^2 _ ^2i„ ^2 ^2_L^2 



^ , , , ^ , , , X 4 mn In rrin — m? In m? m^ + m? „ , , , , 

i^A'(mo,mi) = Fa' (mi, mo) = ;j - — 5 1 - ^i^A(mo,mi) . (15) 

3 L3 rriQ — m| (mg — m|j^ 



For the case with mo ~ mi, we have 

J7 / 2(mo-mi)2 (mo-mi)^ 

FA(mo,mi)« 3 SO^^^ + ' ' ' ' ^^^^ 

FA'(mo,mi)«+ilnm? + ^^^^^ + ... . (17) 
3 bmi 

When all the additional heavy scalar bosons are degenerate m^ = rrijj^ = m^±, we obtain 5$ = 
r$ = 0. In the SM-like limit sin(/3 — q) = 1 with the further assumption tjij^ = m^, we have 

1 m?j-± , , 

In FIG. [l|, we show predictions on the S and T parameters in the THDM together with the 
allowed regions from the precision data for each confidence level. The SM-like Higgs boson mass is 
varied from 117 GeV to 517 GeV (black curve: from up to down), and the SM-like limit sin(/3 — a) = 
1 and m^j- = 300 GeV are taken. We can see that electroweak precision data favor relatively light 
Higgs boson mh < 145 GeV (90% CL). The degenerated mass of the heavy neutral Higgs bosons 
nij^ = is varied from 200 GeV to 400 GeV by the 10 GeV step (dots: from left to right) for 
the given several values of the SM-like Higgs boson mass rrih = 117, 140,240 and 500 GeV. The 
quadratic dependence on the mass difference between additional heavy scalar bosons can be easily 
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FIG. 1: The analysis in the {S,T) plane is shown in the THDM where the SM-like Higgs boson is taken 
to be 117, 140,240 and 500 GeV, with the SM-like limit sin(/3 - a) = 1 and m^± = 300 GeV. The mass 
of heavy neutral Higgs bosons = is varied from 200 GeV to 400 GeV by the 10 GeV step (dots: 



from left to right). Ellipses correspond to electroweak precision limits with 68% (yZSOcr) and 90% {^/A.Gla) 
confidence level. 

understood by the approximate formula for ~ ""^Hi (IH)- Therefore, the deviation of 

the T parameter is insensitive to M. 

In the SM, the mass of the Higgs boson is constrained due to tree level unitarity. It has been 
studied by considering 6x6 scattering matrix of two body scalar states {hh, hz, zz, u^lo^ , hu^ , zoj^) 
where each eigenvalues of scattering matrices are restricted to be less than a criteria ^ as |ao| < 
where is the S wave amplitude matrix. For ^ = 1/2, the Higgs boson mass is bounded to be 



less than about 710 GeV. In the THDM, there are 14 neutral |16], 8 singly charged and a doubly 
charged two body states |1T7|| . In our numerical analysis, absolute values of all eigenvalues for the S 
wave amplitude matrix are required to be less than 1/2 as for a criteria to keep perturbativity[30]. 
For the constraint from vacuum stability, the Higgs potential is required to be positive for a large 
value of the order parameter. In the SM, this condition is expressed by A > at the tree level. In 
the THDM, the condition for vacuum stability is given by [^-20| 



VA1A2 + A3 -hmin [0,A4 - IA5I] > 0. 



(20) 



In FIG. |2|, we show the regions excluded by various theoretical and experimental constraints in 
the THDM parameter space (on the m^-Am plane) assuming the SM-like Higgs boson mass to 
be TUh = 117, 140, 240 and 500 GeV with sin(/3 — a) = 1, where Am = — m^^. The masses of 
neutral scalars and the soft-breaking mass parameter are taken to be degenerate = m'jj = M^. 
Since quartic coupling constants Aj are independent on tan/3 for rrifj^ = M with the SM-like limit, 
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FIG. 2: Theoretical and experimental constraints in the parameter space of the THDM. Uncolored regions 
are allowed by all the constraints we here considered, i.e., tree level unitarity/stability and electroweak 
precision data, and direct search bound of charged Higgs boson, mjj± < 79.3 GeV. The mass and mixing 
parameters are chosen as = rnj^ = m\, with the SM-like limit sin(/3 — a) = 1. In this limit, constraints 
are independent from tan /3. 

unitarity and stability bounds do not depend on tan/3 in the same limit. Regions excluded by 
the conditions from tree level unitarity and vacuum stability are shown as the green and yellow 
areas, respectively, while that excluded by the electroweak precision data at the 90% (68%) CL is 
indicated by the blue (light-blue) area. Although the bounds from perturbative unitarity, vacuum 
stability and the electroweak precision data with the oblique corrections do not depend on the 
type of Yukawa interaction, the direct search results for the charged Higgs boson depends on that 
via the decay process. The region with charged Higgs boson mass below 79.3 GeVpl is shown as 
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the gray area, which is excluded assuming B{H~^ — )• t'^v) + B{H^ — )• cs)=l. Depending on the 
type of Yukawa interaction we may have additional constraints from the flavor physics 

data analyses such as Bs ^,7 (2^, |2|, |2§ , ^ t+v [||] and r lv9{l = e,/i) @ ||]. We 
do not consider these constraints in FIG. |^ because they are model-dependent. In the upper two 
panels for = 117 GeV and m/i = 140 GeV, entire regions of (< 5 TeV) are allowed by all 
the constraints for a relatively small value for |Am|. On the other hand, in the lower two panels 
for m/i = 240 GeV and rrih = 500 GeV, we can see that deviations from the allowed region by 
the electroweak precision data require new contributions to the electroweak precision parameters; 
i.e., a relatively large value of |Am|. Since quartic coupling constants are constrained by the tree 
level unitarity, masses of heavy scalar bosons are essentially determined by the magnitude of 
for Af^ ^ i;^, where the mass difference Am is expressed by Am ~ (Aj — Xj)v'^/M. For a given 
m/i, the magnitude of Am is determined to satisfy the electroweak precision data via the new T 
parameter contribution, which is proportional to (Am)^. Consequently, m^^ is constrained from 
above by unitarity bounds. For rufi = 240 GeV and mh = 500 GeV cases, we have the upper bound 
to be 2 TeV and 1 TeV, respectively. 

We have shown the results in the case with sin(/3 — a) = 1, and rnjj = rnj^ = in FIG. This 
choice of the parameters would be rather special in the sense that there is no tan f3 dependence in 
this case. In Eqs. (|5|) and (P) with sin(/3 — a) = 1, terms dependent on tan/3 are proportional to 
— M^, and then one of Ai and A2 tends to large when tan/3 7^ 1. Therefore, the parameter 
space is more restricted by the unitarity constraints in the case without degeneracy. Also when 
sin(/3 — a) = 1 is slightly relaxed, the bound from tree level unitarity becomes sensitive to tan f3. For 
larger values of tan f3, the bound becomes more restrictive. Consequently, the tree level unitarity 
bound shown in FIG. ^ can be regarded as the most conservative which is independent of the values 
of tan /3. 

Finally, we shortly discuss the implication to the collider phenomenology. At the LHC, a heavy 
SM-like Higgs boson can be found via the gluon fusion gg ^ h 01 vector boson fusion VV — t- h 
{V = W and Z) with the decays into M^VF^*) and ZZ^*^ |32]. The additional neutral scalar bosons 



H and A) would be produced via gluon fusion gg ^ <p [33], associated production with heavy 



quarks pp — )• ttcf}, bbcj) |34], pair production pp — t- W^* — )• (/}H^ [31, B5|, B6l and pp — t- Z* 



AH charged Higgs boson production would be via gb — )• H^t [37], pp — )• W^H^ |38] 

and pp H+H- [pi. As we discussed above, if the SM-like Higgs boson is heavy, a large mass 
splitting between additional heavy Higgs bosons and A (or H) is required. In such a case, 
their decays into a lighter scalar boson associated with a weak gauge boson — )• (j)W^ (or 
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can be significant pO| . These decay modes are kinematically suppressed by the 
degeneracy of scalar boson masses in most of previous discussions, for example, on Higgs boson 



decays in the minimal supersymmetric SM|41|. In addition, the one-loop induced decay process 
— )• W^Z can also be enhanced when the mass difference between A and is large [^]. The 
bosonic decay branching fractions of scalar bosons can dominate over their fermionic decay modes. 
Therefore, detailed studies for these decay modes will be important to test the scenario with large 
mass splitting between additional heavy Higgs bosons fi^ . 

In conclusion, we have analyzed theoretical bounds and experimental constraints in the 
THDM. For a given Higgs boson mass, the magnitude of the mass difference between additional 
heavy scalar bosons can be determined to satisfy the electroweak precision data. However, the 
mass difference requires a large coupling constant in the Higgs potential, and too large coupling 
constant violates tree level unitarity. Therefore, we have found that an upper bound on the 
additional heavy Higgs bosons is obtained when the SM Higgs boson is heavy. For example, 
is bounded to be less than around 2 (1) TeV for m/j = 240 (500) GeV, where = m?^ = m\ 
with the SM-like limit sin(/3 — a) = 1 is taken. Even if the SM-like Higgs boson is found to be 
light (< 140 GeV) our analysis show a possible range of mass splitting in the heavy Higgs bosons 
in the THDM. 
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